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ABSTRACT: The perihelion advance of the orbit of Mercury has long been one of the 
observational cornerstones for testing General Relativity (G.R.). 
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The main goal of this paper is to discuss how, presently, observational and theoretical 
constraints may challenge Einstein's theory of gravitation characterized by (3 = 7 = 1 . To 
£f) ' achieve this purpose, we will first recall the experimental constraints upon the Eddington- 

Robertson parameters 7, f3 and the observational bounds for the perihelion advance of 
Mercury, AiO obs . 

A second point will address the values given, up to now, to the solar quadrupole mo- 
ment by several authors. Then, we will briefly comment why we use a recent theoretical 
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determination of the solar quadrupole moment, J2 = (2.0 i 0.4) 10 7 , which takes 
into account both surfacic and internal differential rotation, in order to compute the solar 
C^) ■ contribution to Mercury's perihelion advance. 

, Further on, combining bounds on 7 and J2 contributions, and taking into account the 

observational data range for Aa> b s , we will be able to give a range of values for (3. 

Alternatively, taking into account the observed value of Acu ;, s , one can deduce a 
' dynamical estimation of J 2 in the setting of G.R. This point is important as it provides a 

solar model independent estimation that can be confronted with other determinations of 
J 2 based upon solar theory and solar observations (oscillation data, oblateness...). 

Finally, a glimpse at future satellite experiments will help us to understand how 
stronger constraints upon the parameter space (7, f3, J2) as well as a separation of the two 
contributions (from the quadrupole moment, Jjj, or purely relativistic, 2oP + 2«7 — (3) 
might be expected in the future. 
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1 INTRODUCTION 



The solar quadrupole moment, J 2 , is one of the fundamental figures in 
solar physics. It provides informations on the distortion of the effective solar 
potential, J 2 being the first perturbation coefficient to a pure spherically 
symmetric gravitational field: 



<&(r,0) = 



00 / F? \ 2 

Ef^) J n P n (cos9) 

n=l ^ ' 



where $ is the solar component of the gravitational potential outside the 
Sun, in polar coordinates (r, 9, 0) with respect to the Sun's rotation axis; 
P n are Legendre functions of degree n. The coefficients J n are thus directly 
related to the distorted shape of the Sun; for instance, for n = 2, J 2 7^ is 
an indicator of the oblatenessQ 

Concerning the Sun, J 2 , which is the most important term, should be used 
as a constraint in the computation of solar models, as the asphericity is a 
probe to test the solar interior. 

Further, detection of long term changes in the solar figure (as there is some 
evidence for J 2 to vary with time) are intended; those have been postulated 
to act as a potential gravitational reservoir that can be a source of solar 
luminosity variations, which in turn, could have significant effects on the 
climate of the Earth ( ||So et al.l9T9[| , |[Roz2001a|| ). 



Today, the quadrupolar moment is also a non negligeable quantity in 
computing the relativistic motion of planets. 

The first time that the solar quadrupole moment was associated with gravi- 
tational motion of Mercury is in 1885, when Newcomb attempted to account 
for the anomalous perihelion advance of Mercury with a modified gravita- 
tional field, manifested by an oblateness Ar [ |Nel895-1898 ]. Indeed, in 1859, 



Le Verrier had observed a deviation of Mercury's orbit from Newtonian's 
predictions, that could not be due to the presence of known planets. But, 
the difference between the equatorial and polar diameters of the Sun of 500 
arc ms, as advocated by Newcomb, was soon ruled out by solar observations. 
And Einstein's new theory of gravitation, General Relativity, could account 
for almost all the observed perihelion advance. 

So, Mercury's perihelion advance readily became one of the cornerstones for 



5 Notice that J2 = — C02, where C02 is the second spherical harmonic coefficient. 
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testing General Relativity; even though, now, a contribution to the perihe- 
lion shift from the solar figure (though less important than first suggested by 
Newcomb) can not be discarded. 

Mercury is the inner most of the four terrestrial planets in the Solar Sys- 
tem, moving with a high velocity in the Sun's gravitational field. Only comets 
and asteroids approach the Sun closer at perihelion. This why the "Mercury 
lab" and minor planets too (see section [4.1. 2| ) offer unique possibilities for 
testing G.R. and exploring the limits of alternative theories of gravitation 
with an interesting accuracy. 

However, the perihelion shift of planets, and hence Mercury, can not be 
measured directly because the perihelion is a Keplerian element whereas the 
motions of the planets are not exactly Keplerian due to mutual gravitational 
interactions and figure effects. So, only an indirect determination can be 
done. One can proceed as follows. The motions of planets, from numerically 
integrated ephemeris, are computed over an interval of time. The time evolu- 
tion of osculating elements is then plot and a polynomial fit of the parameters 
gives the rate of the perihelion advance. If one repeats this procedure in the 
classical Newtonian limit, one gets another set of rates. The difference be- 
tween the two computations, and taking into account the constant general 
precession of the equinoxes, gives the combined effect due to relativistic grav- 
itation and the Sun's quadrupole moment, Au; obs . 

Nevertheless, Au b s depends on how the perturbation elements (for example 
the slow motion of the ecliptic) are taken into account in the computations 



NaRal98o| ; it also depends on the precision and on the data set selected 
from the radar data [|Ral987|l which provide the core of the ephemerides 
computation. 

Furthermore, in the article [Pitl993|| , the author shows that the topographic 
features of Mercury's surface influence the results on the perihelion advance 
inferred from Mercury radar observations. 

These are the main reasons for which the range of Mercury's perihelion ad- 
vance deduced from the radar data remains of a great amplitude (see table 
1, section |7|). 

M. Standish, [|St2000 |, has applied a method analoguous the above mentioned 
method; integrating equations over four centuries, 1800-2200, with and with- 
out the relativistic contribution (the second integration was done by simply 
replacing the speed of light with a very large value). He then computed the 
perihelion of Mercury from both of the runs, one point every 400 days, and 
differentiated the values of the perihelion at each time-point. After fitting 
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a linear function to the differences, the resulting slope from the figure is: 
42.980 ± 0.002 arcsec/cy. Both integrations assumed J 2 = (2.0 ± 0.4) 10~ 7 , 
hence the estimation of the perihelion advance given by M. Standish repre- 
sents solely the purely relativistic contribution. 

At such a level of precision, this solution would scarcely change with future 
ephemeris improvements (or with another set of ephemerides calculations, 
such as those given by the Bureau des Longitudes, since 1889 ||Bureaul989|l ). 

In the following second section, we will describe how, using the most ac- 
curate theoretical value for J 2 , the observed perihelion advance can lead to 
constraints upon the parameters (/3, 7) which describe a generic, metric and 
conservative theory of gravitation. 

In section 3, we will see how, in the setting of General Relativity, Au b s 
can be used to provide a dynamical, solar model independent, estimation of 
J 2- This dynamical value can be confronted with that derived from direct 
measurements of the solar oblateness or indirect ones coming from helioseis- 
mology, which are solar model dependent. 

Finally, we will give, in section 4, an overview of future satellite experiments 
that might be expected to put stronger constraints upon the parameter space 
( (3, 7, J 2) in the future. 

Throughout this article, we will refer to estimated values of Au bs and J2 
from different authors and sources. Those are listed in tables at the end of 
this review, along with the figures, in section [7[ 



2 CONSTRAINTS UPON GRAVITATION 
THEORIES 

2.1 The relativistic advance of the perihelion of Mer- 
cury 

2.1.1 The purely relativistic effect. 

Once correcting for the perturbation due to the general precession of the 
equinoxes (~ 5000 (arcsec/cy)) and for the perturbations due to other planets 
(computed numerically with a Newtonian N-body model: ~ 280 (arcsec/cy) 
from Venus, ~ 150 (arcsec/cy) from Jupiter and ~ 100 (arcsec/cy) from the 
rest), the advance (in regards to the classical Keplerian prediction) of the 
perihelion of Mercury is a combination of a purely relativistic effect and a 
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contribution from the Sun's quadrupole moment. 
It is given by the following general expression?]: 



Au = Au GR 5 
with Au gr = 

S 



(rad/revolution) 

3nR 



a a (1— e 2 ) 

| (2a 2 + 2a 7 - (3) 
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and where the following parameters^ are 
R 



sin 2 i 



1) 



R s 

J2 



the Schwarzschild radius of the Sun, 2G ^ Is 
given in [EPJ2000]; 

the Sun's mass, given in [|A112000|1 ; 



the Sun's radius, given in [EPJ2000 



the quadrupole moment of the Sun for which we take 
the theoretical value of (2.0 ± 0.4) 10~ 7 



(1) 



in 



GodRoz2000l ; 



the semi-major axis of Mercury's orbit, in [[A112000 



the exentricity of Mercury's orbit, in |A112000|| ; 



the inclination of Mercury's orbit, in | A1120Q0 



(2) 

(3) 
(4) 
(5) 

(6) 
(7) 
(8) 



Notice that formula ([!]) is only valid for fully conservative theories. If it is 
not the case, the complete expression is recovered with the following change 



\ (2a 2 + 2a~f - (3) 



R a a(l-e 2 ) 



J 2 (3 sin 2 % — 1J 



+ g (2«i - a 2 + «3 + 2( 2 



" "a 2 " " 



M a M M 
(M a +M M f J 

^3" parametrize preferred- 



, "0J3", a violation of the conservation of the total 



where "Mm" is the mass of Mercury; "«i 
frame effects; and "C2' 
momentum. 

But the extra term is nevertheless negligeable because it is proportional to 
,^ aA ff V 2 ~ ¥¥- ~ 2 10~ 7 (see [|WilU993|| ), and thus negligeable in regards 



In some references, the coefficient of the term containing the contribution of the orbit's 
inclination is improperly written. 

5 Notice that the value of the Schwarzschild radius of the Sun is more accurate than the 
separate values of the gravitational constant, G, and the Sun's mass, M s . 
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to the first one (of the order of unity) or to the second one (of the order of 
10" 4 ). 



"a", "/3" and "7", refer to the Eddington- Robertson parameters of the 
Parametrized Post-Newtonian (P.P.N.) formalism, describing a fully conser- 
vative relativistic theory of gravitation. 

"a" describes the weak equivalence principle; is the amount of non- 
linearity in the superposition law of gravity; and "7" characterizes the amount 
of space curvature produced by unit rest mass. 

The P.P.N, parameters also cover the particular case of Einstein's theory of 
gravitation, General Relativity, characterized bya = /? = 7 = l. 



2.1.2 Constraints upon the Eddington-Robertson parameters. 

First of all, the parameter a is set to unity for any theory that respects 
the weak equivalence principle, well tested (the difference between the accel- 
eration towards the Earth of two test-masses of different composition, relative 
to the sum of those accelerations, is inferior to ~ 10~ 14 . Sec ||Will2001| ). 



Note that the Microscope Mission, selected by the French agency C.N.E.S. 
and scheduled for launch by 2004, has for scientific objective to test the 
equivalence principle up to an accuracy of 10~ 15 , using its well known mani- 
festation, the universality of free fall |[To et al.2000|| . 



Secondly, it is light deflection experiments (measuring the combination 
^p-) that provide so far the best constrains on 7 [|Lebl995|| p|: 

7 = 0.9996 ±0.0017 (9) 



But there is, presently, no independent determination of the parameter 
(3, which appears either in the combination 2a 2 + 2«7 — (3 characterizing the 
perihelion advance, or in the Nordtvedt effect 4(3 — 7 — 3 = rj (see section 

Hi- 



2.2 Theoretical solar quadrupole moment contribution 



6 According to other authors, [Rob et al.1991], the value of this parameter deduced 
from V. L. B.I, measurements is 7 = 1.0002 ±0.00096; whilts Eubanks et al., as quoted by 
|Will200li give i±2 = 0.99992 ±0.00014 (not yet published). 
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2.2.1 The question of the accurate determination of J 2 . 



The evaluation of the solar quadrupole moment, J2, still faces some 
controversy: on one side, the theoretical values strongly depend on the solar 
model used, whereas accurate measurements are very difficult to obtain from 
observations. 

Concerning this last point, let us for example recall some problems: 

(1) the real differences of brightness of the solar limb dependency on the 
latitude; influence of faculae, sunspots and magnetic fields; correlatively, real 
effects due to latitudinal variation of the solar limb darkening function; 

(2) the questioned solar activity (solar cycle) dependency of the Sun's oblate- 
ness. These variations were first conjectured by Dicke et al. in 1985 

IPi ct al. 1987|| . Observations at the Pic du Midi Observatory (France) from 
1993 till 2000 seem to confirm a faint variability reported in previous ob- 
servations made in 1983-1984. Nevertheless, the amplitude of the observed 
variations does not exceed 0.02" — 0.04" over 20 years. (From | Kul998|| . See 
also in ||Roz 1996)1 , fig. 1 and 2; in [[RozR51996|| , and in [[RozRol997|| where the 
authors derive from all the available data a maximum value of J2 of 1 10~ 5 
and an average value of (3.64 ± 2.84) 10~ 6 ); 

(3) and the difficulty to calibrate ground data in regards to atmospheric 
disturbances (local atmospheric refractive indexes and distortions due to at- 
mospheric waves). 

Space experiments have been suggested in order to solve those problems; 
however, first results obtained from the SoHO mission, [|Ku et al.1998 || , have 
established a good concordance with ground-based observations of the oblate- 
ness (and thus J2)- Further comments are found in section iA. 



To illustrate those difficulties, we give a compilation of the main deter- 
minations of J2, based on observations and solar theory, in addition to the 
main critics to the method used (see table 2 and figure 2, section |7|). A 
more detailed historical review can be found in ||Rozl996| or in |[RozRol997| . 
Remark that early estimations of J2, before 1967, using an heliometer of pho- 
tographic plates, often erroneously predicted a prolate Sun (see the second 
table in ||WitDebl987j ). 

In this context, we see that a dynamical determination of J2, using the 
perihelion shift of Mercury, is interesting as it might be confronted to those 
derived from solar model dependent values of the oblateness (see section ^ . 
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2.2.2 The adopted theoretical value of J 2 . 



The theoretical value of J2, used in this article, has been deduced from a 
recent work, where the authors have applied a "differential theory" to a solar 
stratified model, taking into account the latitudinal differential rotation. The 
result is a determination of J 2 as (1.60 ± 0.04) 10~ 7 at the surface of the Sun 
[ GodRozl999al and | GodRoz2000|| ). 



.see 



The value obtained is in agreement with those calculated by Paterno 
Pa et all 996fl , J2 = (2.22 ± 0.1) 10~ 7 , and Pijpers | |Pijl998|| , 



J2 = (2.18 ± 0.06) 10~ 7 , using in their computations the inversion techniques 
applied to helioseismology. 

The slight difference between these values and those of Godier/Rozelot comes 
mainly from the incertitudes on the solar rotation data due to the analytical 
rotation law adopted by |GodRozl999b || (which gives a low velocity rate at 



the equator a bit lower than what is currently observed). But this difference 
does not question the order of magnitude^ 10~ 4 , of the solar contribution in 
Mercury's perihelion advance. This is why we have admitted the theoretical 
range (2.0 ± 0.4) 10~ 7 for J 2 . 

This value can be confronted to the one given by other authors in table 2 or 
figure 2 (section |7|). 

2.2.3 G.R.'s prediction with and without the quadrupole moment 
contribution. 

Using the values of the parameters given in the appropriate references 
(see (§), (|), (|), (§), © and (0)) in (0), plus the value of the period of 
Mercury's orbit given in [|A112000|| , one finds 



6nGM s 

Aloo gr = — ^4 = 42.98 1 arcsec/cy 10 

a 1 — e z )& 



for which the accuracy is on the last digit. 
This is the prediction of the perihelion shift of Mercury in the setting of G.R 
theory, but omitting the contribution of J 2 . 

This raw value is excluded by the last observational data given by |[Anl992 



Excluding the unacceptable estimations, that till recently, reported J2 to be as large 
as 1Q~ 5 (see table 2) , an order o f magnitude larger than the theoretical upper limit allowed 



by lunar librations |EtozBol998 
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SUUUg, [ |Pit200la| , but not by ||Kr et al.l993l and [|Pitl993|| (see table 1, 



section [?D 

But, once the quadrupolar correction is added, using (|8|), this leads to 

Alu gr e [43.000 ; 43.010] (arcsec/cy) for J 2 = (2.0 ± 0.4) 10~ 7 (11) 

which is now consistent with the observations given by ||Anl992|| , ||St2000|[ . 
Pit2001a|| , while still in agreement with [|Kr et al.l993| and [|Pitl993|| (see 



table 1, section |7|). 

This last result also shows that the theoretical prediction for J 2 , argued by 
the authors, is coherent with observations in the setting of G.R. 

An important remark on values adopted for G.R.'s prediction of Mer- 
cury's perihelion advance, Au gr, i n the past is given in |NoWilll986|| . The 
authors also interestingly underline the following fact: "Although of theo- 
retical interest, the difference between these quoted predictions for Mercury's 
perihelion advance has no observational consequence (for present methods of 
evaluation of Mercury's perihelion shift)". Indeed, the predicted general rel- 
ativists contribution to Mercury's perihelion advance is not an input in cur- 
rent procedures testing gravitational theories with the dynamics of Mercury. 
In modern ephemeris used to compute the motion of Mercury, the equations 
of motion already include relativistic post-Newtonian terms which are non- 
periodic. Those contribute to the secular variation of the orbital elements, 
among which, the perihelion of Mercury. The post-Newtonian terms in the 
ephemeris are modulated by a set of parameters (P.P.N, parameters describ- 
ing the gravitational theory, masses or initial conditions of the planets, ...) 
that become part of a multiparameter least-squared fit to the observational 
data (radar, optical data ...) in order to obtain an improved determination 
of the parameters in the least-squares sense. 

However, it is impossible, presently, to fit simultaneously for both the P.P.N, 
parameters and J 2 , the two contributions, relativistic and Newtonian respec- 
tively, to the perihelion shift being too correlated in the case of Mercury 
alone (see section |4.1.2|) . Thus, one can either directly test (fit) the P.P.N, 
parameters assuming a given input value for J 2 in the ephemeris; or assume 
G.R. as the gravitational theory and test J 2 . In the last case, expression ([I]) 
together with ( |TU[ ) are useful to provide a value of J 2 , once a = /5 = 7 = lis 
assumed. 

Nevertheless, the real general relativistic prediction for the perihelion shift 
of Mercury (Au gr) is given unequivocally by (p~0|) , according to present 
values of astrophysical constants. 
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2.2.4 Alternative theories to G.R. gravity. 

General Relativity is often considered today as "THE" relativist theory 
of gravitation. This pure tensor theory corresponds to a Newtonian potential 
that evolves as 1/r, r being the radial coordinate. The theory so far agrees 
with all the observations made in our solar system. Nevertheless, the theory 
of General Relativity can not be the final theory describing gravitation. 

First of all, from a theoretical point of view, General Relativity can not be 
quantified, and this makes it impossible to unify it with other fundamental 
interactions. 

Moreover, the minimal choice of the Hilbert Einstein action, to which G.R. 
corresponds, is not based upon any fundamental principle. Or to express it in 
another way, it is evident that covariance and Newtonian fields approximation 
alone do not determine uniquely the action. Equivalently, nothing guaranties 
that the Newtonian potential is truly universal. 

Any other theory of gravitation would be valid too, as long as it would lead 
to the same predictions as G.R. that have been tested in the solar system, 
with maybe some departures from the Einsteinian theory on larger distance 
scales. 

Also, let us warn that, from the formal point of view, the theory of General 
Relativity is not invariant under conformal transformations. While, if we 
wish to achieve the junction between particle physics, in which conformal 
invariance plays a crucial role, and gravitation, we should consider a theory 
of gravitation that incorporates this property. 

From the experimental point of view, let us notice that General Relativity 
alone still can not reproduce the flat velocity distributions in the vicinity 
of galaxies. The Newtonian potential would indeed predict a decreasing 
distribution. 

We are thus confronted to the following dilemma: either we suppose the 
existence of dark matter, either we modify the potential for galactic distances. 
This second solution would immediately invalidate General Relativity with 
a null cosmological constant. 

Let us also remark that the solution to the "dark matter dilemma" could 
also be a combination of the two solutions cited here above. 

In conclusion, according to the above arguments, it is fundamental to 
conceive that alternative theories to General Relativity, that is to say (5^1 
or/and 7^1, are truly not excluded by the observations... as the case of 
G.R., (5 — 7 = 1, is only a particular spot in the allowed parameter space 
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This is illustrated by plotting ellipses representing the la, 2a and 3a confi- 
dence levels in the 7) plane, owing to Mercury's perihelion advance test, 
for a fixed value of J 2 (See Figure 1 a, b, c). Further, adding constraints on 7 
and (3 coming from tests of the Nordtvedt effect and L.L.R. data (see @ and 
(|T2|)) allows to select a portion of the ellipses in the 7) plane. However, 
G.R always belongs at least to the 3a region in the allowed parameter space 
(P, 7), according to the theoretical bounds on J 2 adopted by the authors (see 
section ggg) . 

Nevertheless, we can conclude that (3 = 1 is not the only allowed case. 

Looking more in details at the contribution of (3 to the perihelion shift, 
we see that the deviation, 1 — j3, from G.R.'s value, is, owing to the error 
bars, of the same order of magnitude as the contribution of J 2 . 
Thus two cases may be envisaged: 

Either (3 < 1, which means that Auj tends to be larger than Aujgr, and the 
effect of 1 — (3 adds to the contribution of J 2 . 

Either /3 > 1, which means that Auj tends to be smaller than Aujgr, and the 
effects of 1 — (3 and J 2 substracts. 

A possible consequence of f3 being different from unity is the Nordtvedt 
effect. Indeed, as soon as the combination of the Eddington- Robertson pa- 
rameters given by 77 = 4/5 — 7 — 3 is non null, the gravitational and iner- 
tial masses of a celestial body are no longer the same (see |[Willl993[1 and 
Will200lj | and references there in) 



New analysis of the Lunar Laser Ranging data (L.L.R.) by [[Willi et al.2001| 
provides rj = +0.0002 ± 0.0009, from which one may deduce the acceptable 
range for j3, using the value of 7 given by (^): 

(3 e [0.9993; 1.0006] . (12) 

This in turn allows us to infer a theoretical shift, Ao;: 

Auj G [42.932 ; 43.057] (arcsec/cy) for J 2 = (2.0 ± 0.4) 10~ 7 . 

We can see that it of course contains the particular case of G.R. (Aujgr), and 
that it is consistent with recent observations (Auj obs ) (see table 1, section |7|). 



Alternatively, owing to the remark made in section |2.2.3| , the fit of the 
most recent ephemeris EPM2000 to accurate ranging observations concern- 
ing the motion of planets (and in particular, the perihelion shift of Mercury), 
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provide an astonishingly precise estimation of the Eddington-Robertson pa- 
rameters (3 and 7 for a given theoretical value of J 2 . 
Indeed, according to reference ||Pit200lE| : 



P = 1.0004 ± 0.0002 and 7 = 1.0001 ± 0.0001. (13) 

for a theoretical value of 2.0 10 -7 for J 2 , in agreement with (j5|). However, the 
uncertainties upon the obtained parameters f3 and 7 are formal deviations, 
and realistic error bounds may be an order of magnitude larger. Moreover, 
this estimation is rather tolerant regarding the assumed value of J 2 . Indeed, 
P = 1.000 ±0.001 and 7 = 1.0005 ± 0.0002 have been obtained using the test 
ephemeris which only differ from EPM2000 by the solar oblateness J 2 = 0.0 
[Pit20Qla||! 



3 Inferring a dynamical value of the solar 
quadrupole moment in the setting of G.R. 

Conversely, one may think to infer the absolute value of the quadrupole 
moment, J 2 , which is necessary (owing the allowed parameter space for (3 
and 7) to be in agreement with observations. But, as mentioned in a remark 
in section 2.2.3, the purely relativistic contribution (2a 2 + 2«7 — (3) and the 



quadrupolar moment of the Sun (J 2 ) are too correlated in the perihelion 
advance of Mercury to lead simultaneously to interesting constraints on (/3, 
7) and J 2 separately. 

This is why, so far, a dynamical estimation of J 2 is made in the setting of 
G.R. 

In the particular case of G. R., the theory parametrized by a — P — 7 = 1, 
we find the results listed in table 3 (section ^) inferred from Au b s (table 1, 
section ^) using equation ([!]). 

Nevertheless, J 2 may not exceed the critical theoretical value of 3.0 10 -6 
according to the argument given in [ IRozBolOOS lI based upon the accu- 



rate knowledge of the Moon's physical librations, for which the L.L.R. data 
reaches accuracies at the milli-arcsecond level. 

Moreover, J 2 has to be positive to be in agrement with an oblate Sun. So, 



8 This estimation does not take into account a possible temporal dependence of J2. 
If such a variability exists, the amplitude is, nevertheless, obviously upper bounded by the 
critical value of 3.0 10 -6 . 
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StOTJ and FltMJIE 



only Standish and Pitjeva's last results 
give interesting dynamical constraints upon J 2 (for the other authors, the er- 
ror bars are too large or J 2 is negative, in contradiction with an oblate Sun). 
Namely: J 2 < 2.89 1(T 7 for the EPM1988 ephemeris model, J 2 < 3.38 1(T 7 
for the DE200 model, J 2 = (1.90 ± 0.16) 10~ 7 for the DE405 model and 
J 2 = (2.453 ± 0.701) 10~ 7 for the more recent EPM2000 model (see table 
3, section |7|). Those values are compatible with the solar model dependent 
theoretical value of J 2 , (|5|), argued by the authors in section [2.2.2| . 



4 Increasing precision in the future 

4.1 Prom Hipparcos to GAIA satellite: towards an as- 
tonishing precision upon (7, (3, J2) 

4.1.1 Light deflection: 7 

Milli-arcsec astrometry is available since 1996 from Hipparcos satellite 
data. The reduction of this data required the inclusion of stellar aberration 
up to terms in v/c 2 , as well as the correction (in ^±2) due to the relativistic 
light deflection in the gravitational field of the Earth and the Sun. Calcu- 
lations for the Hipparcos data were implicitly made in the setting of G.R. 
(a = (3 = 7 = 1), thus allowing for this theory to be checked with a precision 
of 3 10~ 3 on 7. This is of course less accurate than the results based on 
V.L.B.I. measurements |[Rob et a!.1991|| , but Hipparcos opened the door to 



future micro-arcsec astrometry, which can improve the precision upon 7 by 
several orders of magnitude. 

Indeed, in the observational context of light deflection, the satellite GAIA 
GAlA2000"f, one cornerstone of ESA's Space Science Programme, to be 



launched in 2009 (or at least no later than 2012) for a five years mission, will 
increase the domain of observations by two orders of magnitude in length 
(now, light deflection is tested on distances ranging from 10 9 to 10 21 m) and 
six orders of magnitude in mass (now 1 to 10 13 M s ). Moreover, GAIA, im- 
proving Hipparcos' performance, will reduce the avoidance angle towards the 
Sun, thus allowing to measure stronger light deflection effects with a reduced 
parallax correlation. 

This all results into an estimated accuracy of 5 10~ 7 on 7. 

Notice that the quadrupolar moment of the Sun, J 2 , has a contribution 
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to the light deflection that is negligeable in the case of GAIA, owing to its 
non null avoidance angle^. 



4.1.2 Perihelion precession for minor planets: 2a 2 + 2a7 — (3 and J2 

GAIA is also expected to observe and discover several hundred thousand 
minor planets, mostly from the Main Belt. 

All of them will acknowledge a perihelion shift (see equation (|l|)), just like 
Mercury, but with a magnitude in respect to the eccentricity, e, inclination, 
i, and semi-major axis, a, of their own orbit. 

Thus, the relativistic correction per revolution to the orbital motion will only 
be significant for the Apollo, Aten and Armor groups, which means of the 
same order of magnitude^ as for Mercury. (In contrast, it will be about seven 
times smaller for minor planets of the Main Belt). But unlike the Apollo and 
Aten groups, the Armor group are not Earth-Orbit crossers. 
On the other side, the absolute precession rate will be approximately four 
times bigger for Mercury than for members of the Apollo or Aten groups 
owing to their respective revolution periods (and more than 100 times bigger 
for the population of the Main Belt). 

Remark that the perihelion shift of the minor planet Icarus had already been 
used in the past (as early as in 1968 |[LiNul969|] ) in order to infer a dynam- 



ical value for J \. But the non uniform distribution of earlier observations^ 
over the orbit of Icarus and Earth seriously affected the suitability of (just) 
Icarus data in verifying G.R. or estimating J 2 independently (see |Shl965| |, 



Sh ct a!.1968 | and [ |Sh et al.1971 ]). So the estimations of J2 were obtained 



assuming G.R. (see table 4, section |7|). 

The advantages of measuring the perihelion shift of minor planets with 
GAIA, in addition to Mercury's, are multiple. 

First, there will be, of course, an increased precision on individual determi- 
nations of Aa>, due to GAIA's technology but also to the fact that minor 
planets are not as extended as Mercury, and so their position can be mea- 



9 In the case of planets like Saturn or Jupiter being the deflector, the contribution of 
Ji planet to the light deflection effect is non negligeable, due to the important magnitude 
of J2 planet and to the fact that grazing incidence is allowed. 

10 For some exemples sec JGAIA200C| | , page 116 table 1.18. 



11 Based on photographic observations from 1949-1968, for Sh ct al.1971 1, plus 7 
Doppler-shift observations for reference | LiNul969fl ; an d additional observations during 



the encounter with Earth in 1987 for reference 
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sured more precisely. 

Secondly, a statistic on several tens of planets, which is a statistic on Au(a, e) 
(or 5(a, e)), will allow to increase the accuracy on the determination of J2 
and the combination "2a 2 + 2a7 — j3" separately. Remember that those two 
contributions have different dependencies in "a (1 — e 2 )" [|Goul982| . 



Thirdly, by studying the precession of the orbital plane of a minor planet 
about the Sun's polar axis, due to the quadrupolar moment of the Sun (J 2 ) 
but unaffected by relativistic gravitation (2a 2 + 2aj — (3), one should be able 
to dynamically measure J2 independently. This effect being more easily dis- 
cernible for moderately large values of the inclination, i, minor planets like 



Icarus with a large value of i (i ~ 16°) would be truly adequate ( ||Dil965 
||Shl965|| , ||Sh et al. 1968|| ). 

A dedicated simulation still has to be performed to assess the real capa- 
bilities of GAIA in that field. But, so far, an estimation of a precision of 10 -4 
on the combination "2a 2 + 2a7 — /3" from individual determinations of Au, 
seems reasonable; moreover, 10~ 5 should be attainable thanks to statistics 
on several tens of planets. 

From the point of view of J2, GAIA should be more precise than 10 -7 , 
but the accuracy is difficult to assess without an extensive simulation on the 
available sampling of "a (1 — e 2 )". Nevertheless, through measurement of 
perihelion advances, GAIA should provide a more accurate dynamical and 
solar model independent determination of J2, to be confronted with solar 
model dependent predictions from, for example, helioseismology data. 



4.1.3 Resulting constraint on j3 

Using the independent constraint upon 7 obtained by GAIA from light 
deflection, and the constraint upon "2a 2 + 2a7 — (3" from perihelion shifts 
measured by GAIA, one should be able to constraint (3 with a precision of 
3 10~ 4 — 3 10 -5 . This is about two orders of magnitude better than the 
present best determinations due to L.L.R.(see equation (|T^)) or direct fits of 



the data on the detection of 77 [Willi et al.2001 



4.2 Alternative future direct measurements of 7 

GAIA will probably not deliver any results on 7 before the end of its mis- 
sion, but, in the mean time, other space or ground based measurements like 
V.L.B.I.'s, will certainly improve the present determination of 7. See table 5 
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(section |7[) for proposed space missions purely dedicated to the measurement 
of 7. 

As an illustration of further prospects, we can cite the Astrodynamical 
Space Test of Relativity using Optical Devices (ASTROD) Pec et al.2000|| , 
a proposal that has been submitted to ESA in response to a "Call for mis- 
sions proposals for two Flexi-Missions" , but which is not yet accepted. Such 
mission, using time-delay measurements between two spacecrafts orbiting 
the Sun and the Earth, would certainly lead to precisions of the order of 
10~ 6 — 10~ 7 on 7. But if the stability of the clocks/lasers can be reduced 
to KT 18 , then, using the range data of ASTROD as an input for a better 
determination of the solar and planetary parameters, one might dream to get 
a precision of — 10~ 9 on 7! Moreover, from the precise determination of 
orbits, information could be given on the solar quadrupole moment, higher 
moments and j3. Again, providing ultra-stable clocks, precisions of the order 
of 10~ 6 on (3 and 4.5 1CT 8 on J 2 could be reached! 

4.3 A Mercury Orbiter mission: measuring f3 indepen- 
dently from 7 as well as separating (2a 2 + 2cry — m 's 
contribution from 

Scheduled to be launched in 2007 (or 2009) for a 2-5 years mission, Bepi- 
Colombo, has been accepted as E.S.A.'s Cornerstone Mission #5 in 1996 
Bepi2000 l. 

It contains three spacecraft elements, among which a Mercury Planetary Or- 
biter that will considerably help to reduce the error bars on the Eddington- 
Robertson parameters /3, 7 and the Sun's quadrupole moment J 2 . 
Indeed, an estimation of the accuracies attainable has been done thanks to a 
full simulation of radio science experiments with calibration of solar plasma 
noise, non gravitational accelerations and systematic effects. 
The measurement of j3, 7, J 2 and 77 is the output of a complex orbit deter- 
mination process in which radio-metric and calibration acquired during the 
mission are used to provide a complete orbital solution which includes the 
osculating orbital elements of the spacecraft and the planet Mercury, as well 
as the harmonic coefficients of the planet's gravity field (at least to the de- 
gree and order 25). Indeed, precision range and range-rate measurements of 
BepiColombo will constrain the position of the planet's center of mass with 
a precision of about 1 m! Thus allowing a truly precise knowledge of the 
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orbital elements and the secular perihelion shift of Mercury in particular. 
An additional advantage of a Mercury Orbiter, regarding the measure of the 
perihelion shift, is that it would considerably reduce the time scale by com- 
parison to recent data observations that use time averages on a decade time 
scale. This would permit an eventual observation of the variation of the per- 
ihelion rotation due to a possible variation of the solar quadrupole moment 

From the view point of constraining the Eddington-Robertson parameter 
7, time delay measurements of radio signals travelling from the spacecraft to 
the Earth and back, combined with Doppler shift measurements of photons 
should help get a yet more precise determination of 7. BepiColombo should 
then be able to improve over Cassini's expected results (see table 4, section 
0), thanks to frequent solar conjunctions during which the Doppler shift effect 
is maximum. A preliminary analysis of the mission estimates that 7 could 
be accurately determined at 2.5 1CT 6 [Pepi2000|| . 



The precise determination of Mercury's motion would also help measure 
the Nordtvedt effect, rj, with an expected accuracy of 2 10~ 5 . This combined 
with the values found for the perihelion advance (2a 2 + 2a7 — (3) and 7 would 
help lift the degeneracy between (3 and J 2- 

Probing the gravitational field of Mercury at various distances from the 
planet would also help separate the effects of J 2 from those of relativistic 
gravitation (2a 2 + 20:7 — /3), owing to their different dependency in the radial 
distance to Mercury [ |WilU993| . 

Following this idea, advantage can be taken of the large eccentricity of Mer- 
cury's orbit to search for periodic orbital perturbations induced by J 2 and 
relativistic gravity ||Goul982| , ||Willl993|| . 

This is how J 2 would be determined independently by the BepiColombo Mis- 
sion: from the precise determination of the secular nodal precession of the 
planet's orbital plane about the Sun's polar axis, due to the quadrupolar 
moment of the Sun, but unaffected by purely relativistic gravitation. 
All this should lead to a determination of J 2 with a precision of 2 10~ 9 , (see 
lu et aETggg and [ |Bepi200U|| ). 



Notice that the influence of Mercury's topography in determining the per- 
ihelion precession has been stressed by some authors |Pitl993| ], |Pit200lB[ ]. It 
has to be taken into account when processing radar observations, as it might 
help to reduce the systematic errors in ranging. So far, the scarcity of radar 
observational data for Mercury restricts the accuracy of estimates for the 
topographical contribution. But future Mercury orbiters, like BepiColombo 
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or N.A.S.A.'s discovery mission named Messenger^ (to be launched in 2004), 
could remedy to that problem by providing useful complementary data upon 
the topography of the planet. 

Finally, we shall cite an interesting proposal from article |[Tu et al.1996 



(page 24, equation 40), that suggests a measure of (3 independent of 7, testing 
the strong equivalence principle, using a Mercury orbiter on a particular 
resonant orbit. 



4.4 Satellites dedicated to J2 

Future solar probes are expected to determine a more precise value for 

■h- 

Indeed, the quadrupole moment of the Sun can be measured dynamically by 
sending and accurately tracking a probe, equipped with a drag-free guidance 
system, to within a few solar radii of the solar center. J2 is then inferred 
from the precise determination of the trajectory. 

Alternatively, J2 can be inferred from in orbit measurements of solar prop- 
erties. But in this case, the reduction of such a measurement will require a 
better understanding on how solar density models and rotational laws influ- 



ence the multipole expansion of the external gravitational field [|UlHal981 
For example, the micro-satellite Picard is a C.N.E.S.Q mission, due for flight 
by the end of 2005. The expected mission lifetime is 3 to 4 years with a 
possible extension to 6 years. 



The aim of Picard |Pa et al.2001||, is to perform in orbit simultaneous accu- 



rate and absolute measurements of the solar diameter, differential rotation 
and irradiance, in addition to low frequency helioseismology, as a permanent 
viewing of the Sun from a G.T.O. orbit should allow the detection of g-modes. 
Picard should be able to measure J2 with of precision of 10 -8 . 
Notice also that the diameter measurements will be obtained at any latitude 
(sunspots and faculae at limb removed) which should allow the detection of 
a latitudinal variation of the diameter and thus, the quadrupole moment, as 
predicted by the theoretical model used by the authors | |GodKozl999b 



Moreover, Picard will observe the Sun in different wavelength bands, among 



12 But, unlike BeppiColombo, it will not test G.R. nor measure the perihelion pre- 
cession... (see http:\\discovery.nasa. gov\messenger.html and section 9.4 in reference 



Bepi2000 ) 
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Centre National d'Etudes Spatiales (France) 
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which the band used for measurements of the Sun's diameter from the ground, 
535.7 nm. This will permit one to compare space-measurements with ground- 
based ones in order to correct for atmospheric perturbations, and so, to even- 
tually re-calibrate the existing data over the former solar cycle dependencies. 
Thus, the Picard mission is clearly designed to solve some of the problems 
mentioned in section [2.2. 1| for the determination of Ji- 



5 CONCLUSIONS 

We have seen in this article that, so far, the theory of General Rel- 
ativity is not excluded by observations. However G.R. only represents 
one possible point in the still large allowed parameter space 7, J 2 ) and 
alternative theories are also permitted. 
More precisely: 

Future space experiments cited in the last section of this article will con- 
siderably reduce the parameter space in the near future, and so impose an 
even more rigorous test to G.R.. 

As a determination of the solar quadrupole moment is concerned, we 
have stressed the importance of confronting a dynamical deter- 
mination of J2, independent of the solar model (obtained from 
perihelion advances or motion of spacecrafts), to other solar model 
dependent values. 

Presently, this dynamical determination of J2 is still dependent upon the 
gravitational theory (2a 2 + 2«7 — /?), but the future GAIA or Bepi- 
Colombo missions should be able to separate those two effects, 
and so obtain a determination of J 2 that is independent of the 
gravitational theory. 

So far too, the authors can say that their estimated theoretical value of 
J2 = (2.0 ± 0.4) 10~ 7 , which is solar model dependent, together with the 
estimated error bars, is completely coherent with the estimated dynamical 
value resulting from Mercury perihelion advance in the setting of General 
Relativity. 

We conclude by reminding that, if future observations confirm the time 
dependency of J2 (for example, a periodicity with the sunspots cycle), this 
effect will have to be taken into account when using perihelion advance as a 
test for G.R. (or alternatively, to estimate J2). This so far has not been the 
case, as dynamically estimated values of J2 come from a mean over several 



19 



decades (see for example 



Lal992|| ). 
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7 TABLES AND FIGURES 



- Table 1 - 

Inferred correction to G.R.'s prediction for Mercury's perihelion advance 



References 



Perihelion 

Advance 

(arcsec/cy) 



Aii! b s — Auq or 
(arcsec/cy) 



Auj obs /Auj GR 



[|Nel895-1898 



HCU943H 



HCU947H 



Dun 1958] 



[|Wayl966 



3h et al.1972 



MoWarl975 l 



[|5h et al.197^ 1 



HAn et al.1978 



[Brel982a 



[NaRal985] 
[ Brel982b| 



[ Ral987[| 



Kr et al.1986 



EPM1988 
DE200 



[|Ral987|1 



[|An et al.l987H 



|Anl991 j 
[|Anl992 j 



[|Kr et al.l993H 

EPM1988 

DE200 
[Pitl993|1 : 

EPM1988 

DE200 



HSt2Q00 l 
DE405 



Pit2001a 



EPM2000 



~43.37 
42.84±1.01 
42.57±0.96 
43.10±0.44 
43.95±0.41 
43.15±0.30 
41.90±0.50 
43.11±0.21 
43.3±0.2 

45.40±0.05 
45.25±0.05 

42.83 ± 0.12 
42.77 ±0.12 
45.47±0.09 
42.92±0.20 
42.94±0.20 
43.13±0.14 

42.985 ± 0.061 
42.978 ± 0.061 

42.964 ± 0.052 
42.970 ± 0.052 

43.004±0.002 

43.0115±0.0085 



~0.39 
-0.14±1.01 
-0.41±0.96 
+0.12±0.44 
+0.97±0.41 
+0.17±0.30 
-1.08±0.50 
+0.13±0.21 
+0.32±0.2 

+2.42±0.16 
+2.27±0.05 

-0.15 ± 0.12 

-0.21 ± 0.12 

+2.49±0.09 

-0.06±0.20 

-0.04±0.20 

+0.15±0.14 

+0.004 ± 0.061 
-0.003 ± 0.061 

-0.017 ±0.052 
-0.011 ± 0.052 

+0.023±0.002 

+0.0305±0.0085 



~1.01 
0.997±0.023 
0.990±0.022 
1.003±0.010 
1.023±0.010 
1.004±0.007 
0.975±0.016 
1.003±0.005 
1.007±0.005 

1.056±0.001 
1.053±0.002 

0.997 ± 0.003 
0.995 ± 0.003 
1.058±0.002 
0.999±0.005 
0.999±0.005 
1.003±0.003 

1.0001 ± 0.0014 
0.9999 ± 0.0014 

0.9996 ± 0.0012 
0.9997 ± 0.0012 

1.00054±0.00005 

1.00071±0.00020 



Table 1. In the past, planetary motions, necessary to infer the value 
of Au b s , were modeled with classical analytical theories. Presently, 
more precise space experiments require accurate numerical ephemeris. 
Those are made possible thanks to new astrometric methods (radar 
ranging, Lunar Laser Ranging, V.L.B.I.. measurements) that led to 
ranging data uncertainties of only a few meters. 

An interesting review of EPM and DE numerical ephemeris can be 



found in reference [Pit2001b]. It discusses the history of planetary 
motion modeling and describes which data (optical, radar, L.L.R. or 
from space-craft) has been used for the different ephemeris. 

Following this historical path, references [ pNfel895-189"8~l , [ pi!943| ] (cal- 
culated for Julian year J1900.00), [|C11947|1 (J1850.00) and [|Brel982a|1 
(J1900.00) used an analysis of the observed data on Mercury which 
was biased by the assumed theory of gravitation. 
In pi943| J, Punl958j ], |Wayl966 [ and [ |MoWarl975| ], the same old 
analytical perturbation theory as Newcomb's, [Nel895-1898 H, was 



used. Or, in |C1194711 , the results are based on Doolittle's calculations 



of the Newtonian motion, with certain corrections [Dol925|. Those 
methods were not as adequate as present numerical computations. 
While the semi-analytical theory developed in [LesBrel982] did not 
use sufficiently accurate equations of motions of planets. 



See arguments given in [NaRal985] and [Ral987] 



To each reference in table 1 corresponds a value of the advance 
of Mercury's perihelion deduced from observational data: Au! f, s . 
For each of them, one may compute the possible corrective factor, 5, 
to the prediction due to Einstein's Gravity (G.R.), Ac^o GR = 42.981 
arcsec/cy, which does not include the quadrupolar (J2) correction. 
Notice that light deflection measurements constraint 7 to 0.9996 ± 



0.0017 [Lebl995]; while theoretical predictions for the solar quadrupole 
moment, taking into account surfacic and internal differential rota- 
tion, give J2 = (2.0 ± 0.4) 10~ 7 , which means that the solar cor- 



rection to the perihelion advance, 



2.8215 10 



-4 J2 



R a a(l-e 2 ) 



J2 (3 



sin 2 i 



1) 



10- 



is 2.8215 10~ 4 (2.0 ±0.4). 



New L.L.R. data [ [Willi et al.2001| , on their side, provide 
P E [0.9993; 1.0006]. 



Here follow some comments on how the values, Aaj b s , given in table 



22 



1, are inferred from the given references: 

In [p!1943| 1, the discordance O-C between observations and the mod- 
eling theory used by Clemence is — 0".07 ± 0.41. In this theory, 
Clemence took 42". 91 for the general relativistic contribution to the 
precession. Au b s is thus the sum of those two numbers, where 
Clemence's estimation of the probable error in the uncertainty in the 
masses, ±0.6, is taken into account. 

The values given in this table for C11947] result from Clemence's value 
42". 56 ± 0.94 (which is the difference between the total perihelion 
precession observed, 5599". 74 ±0.41, and the Newtonian contribution 
of the other Planets plus the effect of the solar oblateness, 5557". 1 8 ± 
0.85) from which Clemence's erroneous estimation of J2's contribution, 
0".010 ± 0.02, has been removed. 



As reference [ Wayl966 | is concerned, Auj b s is obtained by substract- 
ing the total Newtonian contribution of planets, 531". 26 (that had 
been recalculated by the author of | Way 1966 ] using Marsden's masses 
for planets) and the precession of the equinoxes as calculated by New- 
comb HNel895-1898 f , 5024". 53, from the total precession observed 
5599". 74, as cited from pl943fl . 

The values given for |E>h et al.1972 |] are obtained from their given esti- 
mation of 5 (written in their article as A p with J2 = 0) that includes a 
correction for a typical representation of the topography of the Planet 
Mercury, using equations (|l|) and (|l0|). 

In references [ Brel982"a[ ] and |Brel982b[ |, the perihelion advance had 
been recalculated at J1900.00, using contemporary values for plane- 
tary masses (see [ NaRal985| ] and [Ral987] respectively). For them, the 
constant rate of the perihelion advance due to the equinoxes precession 
is taken to be respectively 5029" .0966 + 2.2223T (cited in pre!982a| ] 
from [[Li et aI.1977f 1) and 5029" .0966 + 2.2274T (cited in ||Brel982b|1 
from [BrcChl981]), where T is in Julian centuries with reference to 
the initial epoch J2000; and the N-body Newtonian contribution from 
planets, 528" .95, from |NaRal985|| . 



Notice that in the following articles, [Kr et al.1986 ], [Kr et al.1993 



[ Pit 1995 1 , the authors improperly writes 42.95 (arcsec/cy) for Au>o gr- 
But the values for A n given in those articles are not affected (see 
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remark in section 2.2. 3j ) . Indeed, the correction to the perihelion mo- 
tion, written as A ir in those articles, has been obtained by fitting the 
EPM88 and DE200 numerical ephemerides to all radar observational 
data (from 1964 to 1989 for Pitjeva's) with the needed parameters (the 
elements of planetary orbits, the radii of planets, the value of AU, etc.). 
The correction to the perihelion motion is thus the "observed" devi- 
ation from the value of Mercury perihelion advanced obtained from 
EPM88 ephemerides with zero J2 assumed on the time interval men- 
tioned. In this present table, we computed Ai0 o b s as Auepmi988+A tt, 
where is Aujepmi988 the mean value of the perihelion advance ob- 
tained from EPM1988 ephemerides, namely, 42". 9806 ~ 42" .981 (arc- 
sec/cy). 

As far as the estimation of the quadrupole moment, J2, is concerned, 
it has been incorrectly done, in those articles, with the value 42". 95 
(arcsec/cy). But it can be recomputed afterwards; assuming G.R. 
(fixed value for a = /3 = 7 = 1) and using the differentiation of for- 
mula ([[]) with the proper value for Aloo GR given in Jig): 

^ 2 = Aoj gr 2-8218 10 ~ 3 ^2 = ^2 PPA/1988 + AJ 2 . 

The rounded up result that appears in those articles are not much 
affected (see our table 3). 



In ||Pit2001b|l , the new estimation of A vr= +0".0055 ± 0".0085 is 
obtained from a fit of EPM2000 numerical ephemeris to radar ob- 
servational data over the period 1961-1997. But now, J2 ppa/2000 = 
2.0 10 -7 was assumed and thus, the mean value of the perihelion shift, 
Aujepm2000 —43". 0060, is different from Au>o GR- 

In [ 5t2000 | , a value of J2 DE405 = 2.0 10 -7 was assumed. But the 
estimated perihelion advance given by M. Standish, 42". 980 ±0".002, 
contains solely the purely relativistic contribution. The quadrupole 
moment's, J2 DP405's, contribution must thus be added. 
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to 






Figure 1. For a given value of J2, the perihelion advance of Mercury constitutes a test of the P.P.N, parameters (3 and 7. 
In the j3 and 7 plane (a set to 1), we have plotted la (the smallest), 2a and 3<T (the largest) confidence level ellipses. 
Those are based on the values for Mercury's observed perihelion advance, Au; a {, s , given in table 1. Notice however that 
the value given by [Nel895-1898 ] as well as those given by [Brel982a], [Brel982b] and [Ral987] have not been taken into 
account. Indeed, the first cited reference did not contain any error bars estimation; the other ones used an improper method 
to evaluate /\& b s (see comments of table 1) and the error bars they provide are truely not realistic ones. Remark also that 
the position of the ellipses varies according to the value of J 2 chosen; but, their orientation is determined by the combination 
(2a 2 + 2«7 — (3) that appears in the expression for /\ui. 

Nevertheless, G.R. is still in the 3<7 contours for the allowed theoretical values of J2 argued by the authors (see section 2.2.2). 
Fig. 1 a, b, c represent the confidence contours for f3 and 7, J2 fixed to its minimum, average and maximum value respectively. 

Additional constraints on f3 and 7 (shaded region) can be taken from the Nordtvedt effect and the L.L.R. data 
(see (12) and (9)). 

They allow to determine a portion of the ellipses which constitute the allowed parameter space for f3 and 7. 
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- Table 2 - 

Estimated values of the Solar quadrupole moment J2 from solar observations and solar modeling 



Year 


References 


Method 


J2 


Critics 


1890-1902 


[AmScl905]* 
[WitDebl987]* 


Direct ground based observation of the solar oblateness 
at Gottingen (heliometer). 


< 4.4 icr 6 


a, c, d 


1891 


[WitDebl987]* 


Rotational theory of the Sun by Tisserand 


< 14 icr 6 




1909 


[WitDebl987]* 


Rotational theory of the Sun by Moulton 


< 20 1(T 6 




1966 


[DiGoll967]* 
[DiGoll974] 
[Di et al. 1986]* 


Direct ground based observation of the solar oblateness 
at Princeton 

(integrated flux from inside till outside the limb). 


(2.37 ±0.23) icr 5 

(2.47 ±0.23) 10- 5 


a, d, e 
f, g, i 




[GolScl968 ] 


Theory of the solar figure obtained from a rotational 
law (based upon stability criteria under differential 
rotation and contemporary surface rotation 
observations) plus a contemporary density model that 
are integrated from the center of the Sun till its 
surface, in order to derive £ at the surface. It is 
further constrained by a solar evolution model. 


< 7.96 lO- 5 


k, 1, n, 

V 


1972 
1973 


[Hi et al.1974]* 
[HiStl975] 


Direct ground based observation of the solar oblateness 
using a F.F.T. edge definition, during periods 
of reduced excess brightness (diameter measurement 
and excess equatorial brightness monitoring). 


(9.72 ± 43.4) 10~ 7 


a, d, f, 
g 


(...) 











Year 


References 


Method 


J2 


Critics 


(...) 










1979 


[Goul982] 


Ground based observation of modes frequency 

splittings in solar oscillation data allow to infer 

an internal radial rotation law from which to deduce J%. 


> 1.2 icr 6 

or 

~ 3.6 10- 6 


a, k, 1, 
p, t, u 
v, x 


1979 


[Hi et al.1982] 


Ground based observation of modes frequency 

splittings in solar oscillation data allow to infer 

an internal radial rotation law from which to deduce J^- 


(5.5 ±1.3) 1(T 6 


a, k, 1, 
p, t, u 

V, w 




[UlHal981] 


Theory of the solar figure obtained from a rotational law 
(based upon a differential rotation model and surface 
rotation observations) plus a density model that are 
integrated from the center of the Sun till its surface, in 
order to derive J2 at thesurface. 


(1.25 ±0.25) 1(T 7 


k, 1, n 




[Kil983] 


Theory of the solar figure obtained from a rotational law 
based upon rigid body like rotation, surface rotation 
observations and a homogenuous density model. 
This provides an upper limit for J2 at the surface. 


< 1.08 10~ 5 




1979 


[Ca et al.1983] 


Ground based observation of modes frequency splittings 
in solar oscillation data allow to infer an internal radial 
rotation law from which to deduce J 2- 


> 1.6 10~ 6 

or 

~ 5.0 10" 6 


a, p, t, 

u, V, X 


1983 


[Di et al.1985] 
[Di et al. 1986]* 
[Di et al.1987]* 


Direct ground based observation of the solar oblateness 
during periods of reduced excess brightness, at Wilson, 
(integrated flux from inside till outside the limb). 


(7.92 ± 0.972) 10" 6 


a, d, f, 


(...) 











Year 


References 


Method 


J2 


Critics 


(...) 










1984 


[Bro et al.1989] 


Ground based observation of p-modes frequency 
splittings in solar oscillation data allow to infer an 
internal angular rotation law from which to deduce Ji- 


(1.7 ± 10%) icr 7 


a, 0, p, 
t, u 


1984 


[Duv et al.1984] 


Ground based observation of p-modes frequency splittings 
in solar oscillation data (...) 


(1.7 ±0.4) 10~ 7 


a, u, v 


1984 


[Di et al. 1986]* 
[Di et al.1987]* 


Direct ground based observation of the solar oblateness 
during periods of reduced excess brightness, 
(integrated flux from inside till outside the limb). 


(-1.53 ±2.36) 10~ 6 


a, d, f, 
g, j 


1985 


[Di et al.1987]* 


Direct ground based observation of the solar oblateness 
during periods of reduced excess brightness, 
(integrated flux from inside till outside the limb). 


(4.72 ± 1.53) HT 6 


a, d, f, 
g, j 


1986 


[Bul986 ] 


Limits on the solar oblateness from the theory of solar 
figure given by Roche's and MacLaurin's models. 
The upper limit of a heavy core is taken to infer J2 ■ 


< 1.1 10~ 5 




1989 


[Dell994] 


Analysis of T. Brown's new helioseismic data 


(7.7 ±2.1) 10~ 6 


a, p, t, 
u 


1990 


[Ma et al.1992]* 


Solar Disk Sextant (S.D.S.): a baloon born experiment 
indirectly measuring the solar angular diameter at 
a variety of orientations using the F.F.T. edge definition 
J2 is then evaluated from the infered solar oblateness 
and from solar surface angular rotation data. 


(±1.68 ±5.70) 10~ e 


b, d, f, 

g, y, Z) 

aa, bb, 
cc, dd, 
ee 


(...) 











Year 


References 


Method 




Critics 


(...) 










1992 


[So et al.1994] 


Solar Disk Sextant (S.D.S.) (...) 


(0.3 ±0.6) 10~ 6 


b, d, f, 

g, y, z, 

aa, dd 


1992-1994 


[El et al.1995] 


Ground based observation of p-mode frequency splittings 
in solar oscillation data obtained from the Birmingham 
Solar Oscillation Network (BiSON) allow to infermm an 
internal rotaion law from which to deduce J 2- 


(2.0 ±0.5) 10~ 7 


a, 1, p, 
t, ii, 
jj 


1992 
1994 


[LySol996] 


Solar Disk Sextant (S.D.S.) (...) 


(1.8 ±5.1) 10~ 7 


b, d, f, 
g, y, z, 
aa, dd 
ff, gg 


1990 
1991 
1992 
1994 


[Pa et al.1996] 


Solar Disk Sextant (S.D.S.) (1992 and 1994) lead to a 
measurement of the oblateness (...) which is used with a 
rotation model in order to evaluate J 2- The surface 
rotation model is constrained by ground based 
observations of acoustic p-modes frequency 
splittings from either the helioseismic network IRIS (1991 
-1992) or BiSON (1992-1994). 


(2.22 ±0.1) 10~ 7 
IRIS 

I 

(2.08 ±0.14) 10~ 7 
BISON 


S.D.S: 

b, d, f, 

g, y, z, 

aa, dd, ff; 
BiSON/ 
IRIS: 
a, p, hh, 
ii 


(...) 











Year 


References 


Method 


J<l 


Critics 


(...) 










1993 
1994 


[R6 et al.1996] 


Direct ground based observation of the solar oblateness 
during periods of reduced excess brightness using the 
distance between both inflexion points of the limb 
profile (scanning heliometer provides diameter 
measurements and excess brightness monitoring). 


(2.57 ±2.36) KT 6 


a, d, f, 


1995 
1996 


[Pijl998] 


Ground based observation of frequency splittings in solar 
oscillation data obtained from the Global Oscillation Network 
Group (GONG) (1995-1996) or space observations of 
oscillations ("a" coefficients) by the Solar Heliospheric 
Observatory (SoHO) (1996) allow to infer an internal angular 
rotation law from which to deduce J 2- 


(2.14 ±0.09) 10~ 7 
GONG 

I 

(2.23 ± 0.09) 10~ 7 
SoHO 

=>■ 

(2.18 ±0.06) 10~ 7 
mean value 


a (GONG), 
k, p 


1996 


[RozR61997] 


Direct ground based observation of the solar oblateness 
during periods of reduced excess brightness using the 
distance between both inflexion points of the limb 
profile (scanning heliometer provides diameter 
measurements and excess brightness monitoring). 


(7.57 ±15) 10~ 7 


a, d, f, 


(...) 











Year 


References 


Method 


J2 


Critics 


(...) 










1996 
1997 


[Ku et al.1998 ]* 


Space observation, by SoHO satellite, of the Sun's full limb 
position (Michelson Doppler Imager -M.D.I.- experiment) 
and brightness allow to infer an oblateness from which to 
deduce J 2 by legendre polynomial fit to the observed limb. 
(During periods of reduced solar magnetic activity). 


(-12.5 ±20.1) icr 7 

I 

(-18.2 ± 17.6) 1(T 7 

=>• 

(-16.8 ± 17.3) 1CT 7 
mean value of 1996-1997 


d, f, q, 
r 




[RozBol998] 


Constraints on J2 from the accurate knowledge of the 
moon's physical librations, for which the L.L.R. data 
reach accuracies at the milli- arcsec level. 


< 3 icr 6 


s 




[Rox2000] 


Theory of the solar figure obtained from a rotational law 
(based upon a differential rotation model -deduced from 
helioseismic inversion- and surface rotation observations) 
plus a density model that are integrated from the center 
of the Sun till its surface, to derive J2 at the surface. 


(2.2125 ±0.0075) lO" 7 


k,l 




[GodRozl999a] 
[GodRozl999b ] 
[GodRoz2000] 


Theory of the solar figure obtained from a rotational law 
(based upon a differential rotation model -deduced from 
helioseismic data and p-modes frequency splittings 
obtained by SoHO- and surface rotation observations) 
plus a density model that are integrated from the center 
of the Sun till its surface, in order to derive J2 at the 
surface. 


(1.6 ±0.4) 10~ 7 


k, 1, m 




[Ku2001] 


Reanalysis of observations. 


2.22 10- 7 





Table 2. To each reference corresponds an estimated value of the solar quadru- 
pole moment, the method used to obtain this estimation (solar observations, solar 
modeling), and some critics we formulate in regards to the method. The year 
given in the table is the date the observations were made (not the date of the 
publication) . 

Notice that some authors^ only provide the value of solar equatorial excess radius 
(Ar = r equ — r po i) in their article. We thus inferred the solar quadrupole mo- 
ment (J2) using the following formula [ Rozl996 | J 2 = 2/3 (Ar — 5r) /vq, where 
Sr = 7.8 ± 2.1 arc ms [RozRol997] is the contribution to J 2 due to the surface 
rotation alone, e = At/tq is the solar oblateness and Tq = 9.6 10 5 arc ms, the 
solar radius (i.e. the best sphere passing through r equ and r po i [ Roz et al.2001 1). 
The critics or remarks made are the following ones: 

(a) Ground based experiments are subject to all kinds of atmospheric perturba- 
tions that have to be modeled. 

(b) Balloon flights are still subject to some differential refraction due to residual 
atmosphere (instability) and problems linked to the stability of the pointing in- 
struments. 

(c) The maximum value for Ar is taken, as the measured minimum leads to the 
erroneous prediction of an oblong Sun. 

(d) Observations of the oblateness have to be done only during periods of reduced 
excess brightness in order to be able to deduce the intrinsic visual oblateness from 
the apparent oblateness obtained with whichever edge definition... until the mech- 
anisms of excess brightness are understood and proper models exists for it. 

(e) Did not take into account the solar surface VT° which could lead to a differ- 
ence in brightness indistinguishable from a geometrical oblateness. 

(f) The choice of the edge of the sun's definition profoundly influences the sensi- 
tivity to excess brightness. 

(The F.F.T. -Finite Fourier Transform- edge definition is highly sensitive to the 
limb darkening shape, but this allow a simultaneous sensitive monitoring of the 
excess brightness, and detecting local/global active regions without reliance on 
solar atmosphere models or other observations.) 

This leads to discrepancies among the different results obtained for oblateness 
measurements made during the same period (even with the same instrument!) 
but using different edge definitions. 

(g) The choice of the edge of the sun's definition profoundly influences the appar- 



14 For those authors mentioned with an asterix in the table, the value of J2 as been 
inferred from their given value of the oblateness, e, or excess equatorial radius, Ar. 
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ent displacement of the Sun's edge attributable to atmospheric seeing. 
(The F.F.T. edge definition is less sensitive to this effect than the Dicke Golden- 
berg integral edge definition.) 

(h) Difficulty to correct for the shift of the inflection point. 

(i) The stated error on Ar for the 1966 and 1983 experiment is a formal stan- 
dard deviation. To make allowance for possible seasonal variations in the locally 
induced atmospheric distortions, the error bars should be increased, possibly to 4 
ms. The 1984 and 1985 results are already corrected for this error and thus the 
derived value of Ji- Notice that 1966 results have often been reinterpreted by the 
authors leading to different conclusions (see [ Dil976| ]). 

(j) Observations in 1983, 1984 and 1985 have been made with a modified in- 
strument (see [ pi et al,1985| ]), by comparison to the 1966 experiment of Dicke- 
Goldenberg, that automatically excluded data that was contaminated by signals 
due to substantial facular patches, as well as color dependent brightness signals. 
The possible existence of a color independent brightness signal is however not 
taken into account. 

(k) Dependent upon the solar density model. 
(1) Dependent upon the solar rotation model. 

(m) Assumes the same rotation rate for the core and for the radiative zone; but 
allows fi(r, 9) to vary with the latitude. 

(n) Uses a model of internal rotation which is assumed to be uniformly differential 
(constant on cylinders) through out the convective envelope, and non differential 
below the convective zone, 
(o) No reliable estimates for the uncertainties. 

(p) Helioseismic data are limited by the error bars to distances above 0.2 R s , near 
the surface and near the poles.. 

(q) A mean limb darkening function is used, while a more realistic model should 
use a more complex function. 

(r) 1996 data set gives noisier results as it was obtained without the active M.D.I, 
image stabilization system. 

(s) Simulations have been performed assuming J2 constant. 

(t) It is difficult to correctly identify individual modes of oscillation for all spatial 

scales. Moreover, the oscillations must be adequately long lived. 

(u) Oscillations, as observed from the ground, can not provide a good measure 

of the rotation rate near the solar pole, because foreshortening limits the viewing 

region. 

(v) Assume that the rotational frequency is independent of the latitude, 
(w) Maximum consistent with the stability of the Sun. 
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(x) Some fits to the data produced rotation curves, fl(r), that were highly un- 
physical. The given value for J2 corresponds to the smoothest curve that fits the 
data. 

(y) A small quantity, the separation between 2 solar limb images, is measured, 
instead of the full solar diameter. This enhances the precision with respect to the 
techniques that measure the full diameter directly, 
(z) The instrument scale can be calibrated for any measurement, 
(aa) The quantity measured is located near the optical axis of the instrument (un- 
like for direct diameter measurements), where the optical system is optimal, 
(bb) No true perpendicular diameters are measured (polar and equatorial). The 
resulting J2 is thus probably less than its real value. 

(cc) The resulting oblateness is ~30° offset from the polar-equator position, 
(dd) Solar photospheric T° and VT° may be a function of the activity solar cycle, 
and so, the use of F.F.T. definitions into data reduction from the S.D.S experiment 
would introduce systematic errors. 

(ee) Gravitational distortions (a non constant wedge angle) of the instrument exist 
that were avoided in the next balloon flights (1992-1994). 

(ff) S.D.S. experiments were made on 2 days, 2 years apart (1992-1994) rather than 
continuously over a period of many years. Moreover, there were no observations 
of solar surface rotation available between 1992-1994. Thus, the large number of 
observations did not allow to lower the uncertainties. 

(gg) Uses a simple model of internal rotation of the Sun as constant angular rota- 
tion on cylinders or on cones. 

(hh) BiSON's helioseismic data imply a solar rotation law which is not compatible 
with that inferred from IRIS's. 

(ii) The rotation model does not take into account helioseismologic observations 
made at different latitudes. 

(jj) It was not possible to find a rotation model that reproduced both the splitting 



data reported by | E1 et al,1995 | and the data from the Big Bear Solar Observatory 
(B.B.S.O.) 

Table 3. This table follows table 1. To each reference corresponds an inferred 
value of the solar quadrupole moment, in the setting of G.R., using the perihelion 
shift of Mercury. 

Notice that concerning references [ Kr et al.1986 |, | Kr et a!.1993f| , [ Pitl993 | and 



Pit2001b], the value of the solar quadrupole moment is calculated according to 



the formula given in the comments of table 1. 
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- Table 3 - 
Inferred solar quadrupole moment 
from the perihelion shift of Mercury, assuming G.R. 



References 



Inferred Quadrupole Moment 



'10- 7 1 



Auj cr 2.8218 10" 



Nel895-1898 



C11943] 



C11947H 



Dun 195^ 1 



Wayl966 



Sh et al.1972 



MoWarl975| 1 



Sh et al.l976[| 



An et al.1978 



Brel982; 



NaRal985|] 



Brel982b|1 



Ral987| 



Kr et al. 19861 
EPM1988 
DE200 



Ral987 



An et al.l987H 



Anl991] 
Anl992j1 



Kr et al. 1993 1 

EPM1988 

DE200 
Pit 1993|| : 

EPM1988 

DE200 



St200Cf| : 



DE405 



Pit2001b| : 
EPM2000 



~ +32.1 
-11.6±83.3 
-33.9±79.2 
+9.8±36.3 
+79.9±33.8 
+13.9±24.7 
-89.1±41.2 
+10.6±17.3 
+26.3±16.5 

+199.4±4.1 
+187.1±4.1 

-12.3 ± 10.0 
-17.1 ± 9.7 
+205.2±7.4 
-5.0±16.5 
-3.4±16.5 
+12.3±11.5 

+0.33 ± 5.03 
-0.25 ± 5.03 

-1.40 ± 4.29 
-0.91 ± 4.29 

+1.90 ± 0.16 

+2.453 ± 0.701 
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- Table 4 - 
Inferred solar quadrupole moment 
from the perihelion shift of Icarus, assuming G.R. 



References 



LiNul969 


Lal992 


1 



Inferred Quadrupole Moment 
h 

(by fitting the parameters) 



(+1.8 ±2.0) 1(T 5 
(-0.65 ±5.84) 10- 6 
or < 2 10~ 5 



Table 4. To each reference corresponds an inferred value of the solar quadrupole 
moment, in the setting of G.R., using the perihelion shift of Icarus. 
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h 

0.00008 - 
0.00006 - 
0.00004 - 
0.00002 - 



1900 1920 1940 + 1360 



I 13( 




★ Sun 



a Mercury 



Icarus 



Years 



Figure 2. Different estimated values of the solar quadrupole moment, J%, versus the date when the respective observations 
were made. There are 3 types data points: values estimated from solar models and observations (table 2), values infered 
from the perihelion shift of Mercury (table 3) and those obtained from Icarus' (table 4). 



- Table 5 - 
Proposed space missions dedicated to 7 



References 


Method 


Mission 


Expected 
precision 
on 7 


[BeGil998 ] 


doppler measurement of 


Cassini launched in 1997, 


10" 4 - 10" 5 


[Ie et al.1999] 


the Solar gravitational deflection, 
the first time this method is used 


experiment in 2002-2003 




[Fi et al.1995], [GAIA2000] 


relativity gyroscope experiment, 


Gravity Probe B (2002) 


6 10" 5 


http://einstein.standford.edu 


geodetic precession measurement 






[Bepi2000], [Tu et al.1996] 


output of orbit determination, 


Mercury Orbiter, within 


2.5 10" 6 


http:/ /www.estec. esa.nl/ 


time delay and doppler shift 


BepiColombo (2007/2009) 




spdwww /future /html / 


measurements 






meo2.htm 


(see section 4.3) 






[Ref999] 


Projet d'Hologe Atomique par re- 


International Space Station 


1 1(T 5 


http:/ /www. cnes.fr /WEB 


froidissement d'Atomes en Orbite 


(IIS) (2004/2005) 




_UK/activities/index.htm, 


(PHARAO clock), a swiss hydro- 






see in "Understanding the 


gen maser clock to provide a long 






universe" , "Fund. Phys." 


term frequency standard, associated 
with the IIS, to form the Atomic 
Clock Ensemble in Space (ACES) 






[Rel999] 


Time delay / light deflection measu- 
rements (see section 4.3) 


Solar Orbit Relativity Test 
(SORT) (after 2010) 


1 10 -7 



Table 5. Direct measurement of 7: Some space experiments dedicated to 7. 
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